Astronomy & Astrophysics manuscript no. aa3546revc 


February 5, 2008 


(DOI: will be inserted by hand later) 





An empirical temperature calibration for the Aa photometric 

system. I. The B-type stars* 

E. Paunzen 1 , A. Schnell 1 , H.M. Maitzen 1 



in 
o 
o 

(N 

Or 
<D 

(N 



> 
On 

o 

On 

o 
in 
o 

S3 

a,: 

6 

CO 



X 



Institut fur Astronomie der Universitat Wien, Tiirkenschanzstr. 17, A-1180 Wien, Austria 
Received 1 June 2005 / Accepted 21 August 2005 

Abstract. We establish an empirical effective temperature calibration of main sequence, luminosity class V to III 
B-type stars for the Aa photometric system which was originally developed to detect magnetic chemically peculiar 
objects of the upper main sequence (early B-type to early F-type) at 5200 A. However, this system provides the 
index (gi — y) which shows an excellent correlation with (B — V) as well as (6 — y) and can be used as an 
indicator of the effective temperature. This is supplemented by a very accurate color-magnitude diagram, y or V 
versus (gi — y), which can be used, for example, to determine the reddening, distance and age of an open cluster. 
This makes the Aa photometric system an excellent tool to investigate the Hertzsprung-Russell-Diagram (HRD) 
in more detail. Using the reddening-free parameters and already established calibrations within the Stromgren 
uvbyP, Geneva 7-color and Johnson UBV systems, a polynomial fit of third degree for the averaged effective 
temperatures to the individual (gi — y)o values was derived. For this purpose, data from the literature as well as 
new observations were taken resulting in 225 suitable bright normal B-type objects. The statistical mean of the 
error for this sample is 238 K which is sufficient to investigate the HRD of distant galactic open clusters as well 
as extragalactic aggregates in the future. 
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1. Introduction 

One of the most important observational diagnostic 
tools of astrophysics is the Hertzsprung-Russell-Diagram 
(HRD) allowing the study of the correlation between the 
effective temperature and the absolute magnitude (or lu- 
minosity) of astronomical objects. Virtually all stellar as- 
trophysical models are tested according to the HRD. 

The absolute magnitudes of stellar objects can be 
derived directly via parallax measurements, appropriate 
photometric indices (e.g. Stromgren j3) or on a statisti- 
cal basis in open and globular clusters. The errors of such 
combined estimates are already rather small (< 0.1 mag). 

The photometric calibration of effective temperatures 
is still a very tricky business with several pit falls. The 
applied calibrations depend on the investigated spectral 
range and the physics introduced in the models. Smalley 
& Kupka (1997) give an excellent overview of how con- 
vection, for example, can lead to significant deviations for 
A-type and cooler stars. The statistical calibration of effec- 
tive temperatures via photometric indices has been done 
since the introduction of photometric systems (Johnson 
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1958, Stromgren 1966 and references therein). The ap- 
plied calibrations have become more precise and sophisti- 
cated as the theoretical stellar atmospheres and the input 
physics became more realistic. Furthermore, the amount 
of available photometric and spectroscopic data is increas- 
ing. 

Almost thirty years ago, Maitzen (1976) introduced 
the narrow-band, three filter Aa photometric system in 
order to investigate the flux depression at 5200A which is 
very likely enhanced by the effects of magnetic radiative 
transport phenomena within magnetic chemically pecu- 
liar stars of the upper main sequence (Kupka et al. 2003). 
Since then, more than two dozen papers have presented 
photoelectric and CCD results for galactic field stars, open 
clusters as well as the Large Magellanic Cloud in the rele- 
vant spectral range from early B to early F-type (Maitzen 
1993, Paunzen et al. 2003, Paunzen et al. 2005a,b). The 
Aa system is based on the three narrow band filters g\ , gi 
and y from which two indices are calculated (g\ — y) and 
a = gi — (<7i +y)/2. An a versus (g\ — y) diagram can eas- 
ily detect magnetic chemically peculiar stars (Sect. 12.111 
whereas the (gi — y) index shows an excellent correlation 
with (B — V) as well as (b — y) and can be used as an 
indicator of the effective temperature. The main result 
of detecting peculiar stars is supplemented by a very ac- 
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curate color-magnitude diagram, y or V versus (g\ — y), 
which can be used, for example, to determine the red- 
dening, distance and age of an open cluster (Claret et al. 
2003). 

In the last few years, much theoretical effort was made 
to explain the results of the Aa photometric system in 
astrophysical terms. Kupka et al. (2003, 2004) established 
a synthetic photometric Aa system and confirm the ob- 
served dependency of the a index as a function of var- 
ious colour indices sensitive to the effective temperature 
and surface gravity variations within the Stromgren uvbyfi 
and Johnson U BV photometric systems using fluxes from 
ATLAS9 model atmospheres as well as most recent atomic 
line data together with opacity distribution function for 
individual chemical compositions. Furthermore, Claret et 
al. (2003) calculated isochrones, taking into account mass 
loss during the main sequence evolution, for the Aa photo- 
metric system on the basis of modern equations of state in- 
cluding partial ionization through Saha's approximation, 
the pressure of gas and radiation as well as the equations 
for degenerate electrons. 

We now propose to establish an empirical effective 
temperature calibration for main sequence (luminosity 
class V to III) stars in terms of (g± — y)o. This is most 
important in studying very distant galactic open clusters 
and extragalactic systems for which, in general, no photo- 
metric data within a standard system are available. The 
absolute magnitudes and thus luminosities can be easily 
estimated via y and the appropriate isochrones. In this 
first of two papers we deal with B-type stars because : I) 
There are reddening-free parameters available allowing us 
to check the dereddening procedure for (gi — y) and 2) 
The stellar interiors for these objects are very similar over 
the whole spectral range. The second paper will deal with 
A-type to early F-type objects for which the calibration 
has to be more sophisticated due to the increase of line 
blanketing and luminosity effects. 

For the empirical temperature calibration, we use a ho- 
mogeneous sample of bright (V < 7 mag), apparently nor- 
mal type objects from the literature. But we also present 
new observations for 99 B-type objects that are included 
in our analysis. In total, 225 stars were used to derive ef- 
fective temperatures within the Stromgren uvbyfl, Geneva 
7-color and Johnson UBV systems which were then ap- 
plied to establish a calibration in terms of (gi — y) . The 
final calibration is valid for effective temperatures between 
33000 and 10000 K and yields a statistical mean error of 
238 K for the whole spectral range. 

2. Sample of program stars 

We have selected "normal" B-type objects fulfilling the 
following criteria: 

— classification as B-type, luminosity class V to III 

— no significantly deviating Aa values, i.e. ±I0mmag 

— not listed in the catalogue of Ap and Am stars by 
Renson et al. (1991) 
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Fig. 1. The distribution of Johnson V and (B — V)q, see 
Sect. |3| for our sample of 225 main sequence B-type ob- 
jects. 



— available data within the Johnson, Stromgren, Geneva 
and Aa photometric system 

Binary systems of all kinds and high v sin i stars 
were a-priori not excluded. The Johnson, Geneva and 
Stromgren colors were taken from the General Catalogue 
of Photometric Data (GCPD, Mermilliod et al. 1997). The 
Aa photometry is from Vogt et al. (1998) with additional 
measurements as described in Section 12.11 The (gi — y) 
data published by Vogt et al. (1998; see Table 4 therein) 
are not available in electronic form but were kindly pro- 
vided by the authors. 

The following stars have inconsistent photometric mea- 
surements and were therefore excluded: HR 345, HR 1375, 
HR 1617, HR 2870, HR 3470, HR 8854 and HR 8887. 
From the available data, we are not able to decide if mea- 
surement errors, a wrong identification or a binary nature 
causes these discrepancies. 

The final list comprises 225 objects that satisfy 
our criteria. The complete list of program stars is 
only available in electronic form at the CDS via 
anonymous ftp to cdsarc.u-strasbg.fr (1 30.79.1 25.5), 



http://cdsweb.u-strasbg.fr/Abstract.html or upon request 
from the first author. This table includes the identification 
of objects, the (gi — y)o, (B — V)o, X and (u — b) values, 
V magnitudes, effective temperature with the correspond- 
ing errors (Fig. 0] and Sect.0J, wsini values and spectral 
types, respectively. 
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The distribution of Johnson V and (B — V)q for our 
sample is shown in Fig. ^ (see Sect. [21 for the estimation 
of the reddening). The peak of the (B — V)q values is at 
about —0.10 mag (B8, Table ^) which reflects the coinci- 
dence that the chemically peculiar stars which were mea- 
sured within the Aa photometric system also peak at this 
effective temperature (Schneider 1993) and the normal- 
type objects, used in this investigation, served as standard 
stars. 

2.1. Additional observations 

We observed a sub-sample of the bright B-type objects 
listed in Cowley (1972) in the Aa photometric system. 
She determined MK spectral types for all stars classified 
as B8 in the Bright Star Catalogue north of —20° in order 
to investigate the appearance of CP objects. The results 
of the Aa photometry are discussed in Sect. 12.21 

These photoelectric measurements in the Aa photo- 
metric system were carried out by one of us (A. Schnell) 
at the 60 cm telescope of the Leopold- Figl-Observatory on 
Mitter-Schopfl (880 m above sea level, 40 km southwest of 
Vienna) . The telescope was equipped with a single channel 
photometer and a thermoelectrically cooled EMI 9844A 
photomultiplier using a diaphragm of 35" (only exception- 
ally 47"). The integration times per filter ranged from 15 
to 25 seconds. The g\ and gi filters are identical to those 
described in Maitzen & Floquet (1981) while the third 
one is a conventional Stromgren y filter (Crawford 1978). 
The observations were obtained over 26 nights between 
October 1991 and August 1996. 

The photometric measurements were reduced in the 
usual way, correcting for airmass and extinction. The pro- 
gram stars were observed continuously over the years (up 
to 35 times) which allows us to correct for different zero 
points and instrumental effects. The errors of the means 
for a and (gi — y) are in the range between 1 and 7 mmag, 
with a "mean error" of 3 mmag only. This proves the high 
quality of the individual measurements and supports the 
reduction process. 

Three stars (HR 1363, HR 1617 and HR 8821) are 
in common with the paper by Vogt et al. (1998). The 
measurements of (gi — y) and a agree within 2 mmag for 
these objects. 



2.2. Aa photometry 

The Aa intermediate band photometric system samples 
the depth of the 5200A flux depression by comparing the 
flux at the center (A c = 5210A, FWHM= 120 A, g 2 ), with 
the adjacent regions (5020A, 120A, g x and 5500A, 230A, 
y). The respective index was introduced as: 



a = .92 - (gi + y)/2 



(1) 



It is optimized to detect magnetic chemically peculiar 
(CP2 and CP4) stars, but is also capable of detecting a 
certain percentage of non-magnetic CP1 and CP3 objects 
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Fig. 2. Observed a versus (6 — y)o diagram for our sam- 
ple of B-type objects. The solid line is the normality line 
whereas the dotted lines are the confidence intervals cor- 
responding to 99.9%. 



(Paunzen et al. 2005a). The nomenclature of the differ- 
ent CP groups is according to Preston (1974). The group 
of classical Be/shell and metal- weak stars (Paunzen et al. 
2002) can be investigated. Since a is dependant on the 
temperature, the intrinsic peculiarity index had to be de- 
fined as 



Aa = a-a Q [(b~y);(B~V);(g 1 -y)} 



(2) 



i.e. the difference between the individual a- values and 
those of non-peculiar stars for the same color. The locus 
of the ao-values has been called the normality line. 

Figure [3 shows the results of our observations We in- 
dicate the different CP groups and known Be/shell stars 
with individual symbols. The normality line was calcu- 
lated using (b — y)o of all "normal" type objects, i.e. 
those not included in the catalogue of peculiar objects 
by Renson et al. (1991): 



a = -0.120(1) + 0.130(9) • (b - y) 



(3) 



Known binary system were a-priori not excluded. The 3 a 
limit around the normality line is between 3 and 6 mmag. 

As expected (Paunzen et al. 2005a), the magnetic CP 
stars show the most significant positive Aa values. In 
general, CP3 stars tend to lie above the normality line 



4 



E. Paunzen, A. Schnell, H.M. Maitzen: An empirical temperature calibration for the Aa: B-type stars 



whereas Be/shell stars are below it. There are two ex- 
ceptions, HR 6664 and HR 7721 which might be in their 
shell phases for which positive Aa values are expected 
(Maitzen & Pavlovski 1987). We note that the two ap- 
parently magnetic CP stars HR 481 (B8 IIIp Si) and HR 
7401 (B8 IV He- weak) show no significant positive Aa 
values and should therefore be considered as misclassified. 
Only one apparently normal type object, HR 7737 (B9 
IV-V, Aa = +14mmmag), deviates significantly from the 
normality line. This star is a known spectroscopic binary 
system with unresolved components of 7.3 and 7.5 mag 
(combined magnitude: 6.79 mag; Germain et al. 1999). 
This seems to influence the photometric values in the sense 
that it mimics a peculiar object. 

3. The estimation of the reddening 

The reddening for B-type stars within the solar neighbor- 
hood is, in general, estimated using photometric calibra- 
tions in the Stromgren uvbyfi (Crawford 1978) and the 
Q parameter within the Johnson UBV system (Johnson 
1958). These methods are only based on photometric in- 
dices and do not take into account any distance estimates 
via parallax measurements. 

The reddening in the Stromgren uvby/3 photomet- 
ric system is based on the comparison of the reddened 
(b — y) and c\ with the unreddened (it — b) and (5 in- 
dices (Crawford 1978). The procedure of the Q method 
is straightforward and has been described in detail by 
Johnson (1958). Here are the basic correlations from this 
reference: 

Q = (U-B)- 0.72 ■ (B — V) — 

Om-(B-V) 2 (4) 
E{B -V) = (B — V) — 0.332 • Q (5) 

Throughout this paper, we use the following relation: 
Ay = 3AE(B -V) = A.3E(b-y). Figured] shows the com- 
parison of the derived Ay values for both methods. The 
agreement is very good indicating that photometric mea- 
surements within the Johnson and Stromgren systems are 
consistent. Otherwise, a severe deviation from the linear 
correlation would occur. The distribution of the adopted 
values (mean of both methods) shows that most of the 
program stars have an absorption which is below 0.1 mag. 

However, we also derived the interstellar reddening us- 
ing the model proposed by Chen et al. (1998) who com- 
bined galactic reddening maps, which are derived from 
open clusters as well as from galactic field stars with pub- 
lished empirical reddening laws from the literature. As 
input parameters, the galactic coordinates and the dis- 
tance from the Sun (i.e. derived from Hipparcos parallax 
measurements) are needed. The error of the latter severely 
influences the error of the derived reddening. For program 
stars with distance errors smaller than 15%, a very good 
agreement with the results of the photometric calibrations 
has been found. 




0.0 0.2 0.4 0.6 0.8 1.0 
A y (uvbyP) [mag] 




0.0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 
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Fig. 3. The correlation of the absorption Ay estimated 
from the Stromgren uvby(3 and the Q method (upper 
panel) . The statistical error for both methods is indicated 
with the symbol in the lower left corner. The distribution 
of the adopted Ay values is shown in the lower panel. 



4. The calibration of the effective temperature 

To reach our final goal, an empirical effective tempera- 
ture calibration of (gl — y)o for B-type stars, we calibrated 
this astrophysical parameter for our sample using the pub- 
lished calibrations in the Geneva, Stromgren and Johnson 
photometric systems. Those calibrations are derived inde- 
pendently of each other which allows us to detect possible 
inconsistencies due to, for example, spectroscopic binary 
nature. We will now discuss the applied calibrations in 
more detail. 

Geneva system: detailed calibrations were published by 
Cramer (1984, 1999) and Kiinzli ct al (1997). They are all 
based on the reddening-free parameters X and Y which 
are valid for spectral types hotter than approximately AO. 
The results are therefore independent of the estimation of 
Ay for our program stars. 

Stromgren system: the most recent and widely used ref- 
erence is Napiwotzki et al. (1993). For stars hotter than 
11000 K, the unreddened [u — b] and for cooler B-type ob- 
jects, the ao index are used to calibrate the effective tem- 
perature. The latter is not reddening free. 
Johnson system: we calculated the Q values for luminos- 
ity class III and V objects according to the Tables listed 
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Fig. 4. Mean relation between the effective temperature 
and (u — b), X, (B — V)q as well as (gi — y)o for B-type, 
luminosity class V to III objects. 



by Schmidt-Kaler (1982). The (B - V) values for those 
luminosity classes were transformed into effective temper- 
atures using the results by Code et al. (1976, Table 7). 
As final correlation we derived a polynomial fit for depen- 
dence of the effective temperature on Q as 

logT e// = 3.983(5) - 0.31(3) • Q + 0.33(3) • Q 2 (6) 

which is valid for all B-type, luminosity class III to V ob- 
jects. This relation has to be treated as an averaged statis- 
tical result. If one uses the individual measurements from 
Code et al. (1976) for the ten stars that are in the relevant 
effective temperature and luminosity regime (Sect.[2J), the 
following relation is derived: 

logT e// = 3.994(29) - 0.25(15) • Q + 0.38(15) ■ Q 2 (7) 

A comparison of these two relations give 
AT eff = +141(147) K. 

The individual effective temperature values within 
each photometric system were first checked for their in- 
trinsic consistency and then averaged. These final values 
together with the standard deviations of the means are 
listed in the electronically available table. No outliers in 
any photometric system were detected. 

5. The calibration procedure for the Aa 
photometric system 

The starting point is one set of (gi — y) measurements, for 
example, of an open cluster. The zero points of data sets 
from different instruments and filter systems may vary. 
An overview of all currently used Aa filters is listed in 



Paunzen et al. (2005a, Table 2). They have also inves- 
tigated the effects on the absolute photometric values, 
which arc in the range of the observational errors for the 
photoelectric as well as CCD technique (a few mmags 
only). 

The reddening coefficient k of the relation E(gi — 
y) = k ■ E(B — V) was estimated which can then be trans- 
formed to other photometric systems. We derived the cor- 
relation coefficients for (B—V) versus (gi—y) and then ap- 
plied these coefficients to (B— V)o using the estimated red- 
dening for each individual star as well as {g\ —y) to derive 
(gi—y)o and thus E(gi—y). Our final value, k — 0.39(2) is 
in excellent agreement with the values derived from open 
cluster data (see Maitzen 1993 and references therein). 

We then dereddened the individual indices and estab- 
lished the linear correlations between those parameters. 
We define a "standard (<?i — y) system" which is set to 
(B - V) = {b- y) = {gi - y)o = as follows: 

(B - V) = -1.413(7) + 2.82(1) ■ ( 9l - y) (8) 
(b - y) = -0.664(5) + 1.333(9) • ( 9l - y) (9) 

This system will also be kept for the second paper. 

Figure 0] shows the relation between the effective tem- 
perature and the different temperature sensitive indices 
for the four investigated photometric systems. The tem- 
peratures range from approximately 33000 to 10000 K 
which covers the main sequence B-type stars. However, 
there are only four stars with temperatures hotter than 
25000 K (B1.5V). We checked the temperature relations 
with and without those four data points and found no 
significant differences. Our final calibrations are: 

logT e// = +4.520(3) - 0.754(15) ■ (u - b) + 

+0.413(2) • (it - b) 2 - 0.108(11) • {u - 6) 3 (10) 
= +4.546(3) - 0.698(13) • X + 

+0.367(17) -X 2 - 0.092(6) ■ X 3 (11) 
= +3.956(5) - 1.04(7) • (B - V) + 

+2.89(24) • (B - V)l (12) 
= +3.909(7) - 6.47(48) • {g x - y) + 

-47(9) ■ ( 9l - y)\ - 425(60) • ( 9l - yf (13) 

with the mean of the errors for the whole sample of 
AT e/ f [(« — b), X, (B — V) , (51 - y)o]= [157,146,333,238K] . 
These are statistical errors and should be treated as such. 
The mean relations between the effective temperature and 
(B — V)o, (u — b), X as well as (gi — y)o depending on 
spectral types are listed in Table The agreement with 
standard values for those indices from the literature (Code 
et al. 1976, Crawford 1978, Cramer 1999) is excellent. 

It is known that high rotational velocities can alter 
the photometric indices significantly (Collins et al. 1991). 
The break-up velocity ranges from about 540 to 400 kms -1 
for B0 to B9 stars (Townsend et al. 2004). However, the 
inclination i is the crucial point for the comparison of 
these models with observations. The prototype AO star 
Vega has a very low vsini of 22kms~ 1 but an equatorial 
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100 200 300 400 100 200 300 400 



v sim [kins' 1 ] v sim [kms 1 ] 

Fig. 5. No systematic correlation of AT e // and the vsmi 
values were found for (u — b), X, (B — V)q as well as 
(.91 - V)a- 



velocity of 160 kms -1 (Hill et al. 2004). We searched in 
the papers by Glebocki & Stawikowski (2000) and Abt 
et al. (2002) for available wsini values for our program 
stars, and averaged individual measurements from the dif- 
ferent sources. In total, there are measurements for 175 
stars of our sample available with the highest value of 
390 kms -1 for HR 3502. Figure shows the diagrams of 
AT e /j =T e f f(orig) —T e ff(calib) versus i>sini for the dif- 
ferent indices. There is no obvious correlation between 
these two parameters in any of the diagrams. The effect 
of the rotational velocity is therefore less significant than 
the precision of the method itself. Thus, the well-known 
alteration of stellar colors caused by high rotational ve- 
locities (Collins et al. 1991) cannot be distinguished from 
the overall statistical errors resulting from the calibration 
process. 

The overall procedure of deriving the effective temper- 
atures for main sequence B-type objects in the Aa photo- 
metric system should be: 

— Estimate the reddening, for example via isochrones for 
open cluster members 

— E( 9l -y) =0.39 -E(B-V) 

— Apply the reddening correction for all individual in- 
dices 

— Transform the (g\ — y)o via the standard relations 

— Check the intrinsic consistency of all available mea- 
surements according to the spectral type - effective 
temperature - photometric indices correlation 

— Apply the effective temperature calibration 

The estimated effective temperature from the Aa system 
can in addition be compared with calibrated values from 
other photometric indices. 



Table 1. Mean relation between the effective temperature 
and (B — V)o, (u — b), X as well as (gi — y)o for B-type, 
luminosity class V to III objects. 
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6. Conclusion and outlook 

An empirical effective temperature calibration for the Aa 
photometric system for main sequence (luminosity class 
V to III) B-type stars was established. This system orig- 
inally was developed to search for magnetic chemically 
peculiar objects of the upper main sequence. However, it 
also provides the temperature sensitive (gi — y) index as 
well as y which can be used to determine the temper- 
ature and absolute magnitude (or luminosity) of a star. 
For this purpose, data from the literature as well as new 
observations were taken resulting in 225 suitably bright 
normal type objects. Known binary stars and high wsinj 
stars were a-priori not excluded. All program stars were 
first dcrcddcncd and then calibrated using already estab- 
lished methods in the Stromgren uvby(3, Geneva 7-color 
and Johnson UBV photometric systems. These calibra- 
tions are based on (u — b), X and (B — V)o, respectively. 

The intrinsic consistency for the individual stars was 
checked and as final step, the averaged effective temper- 
atures were calibrated in terms of (gi — y)o yielding a 
polynomial fit of third degree. The statistical mean of the 
error for the sample of 225 stars is 238 K which is very sat- 
isfactory and sufficient to investigate the HRD of distant 
galactic open clusters as well as extragalactic aggregates. 

The second part of this series will investigate an em- 
pirical effective temperature calibration of A-type to early 
F-type objects for which the calibration has to be more 
sophisticated due to the increase of line blanketing and 
luminosity effects without the availability of an indepen- 
dent reddening free parameter for the Johnson UBV and 
Geneva 7-color system. 
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